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'H: ABSTRACT 

We consider an addition of the term which is a square of the scalar curvature to 
the Einstein-Hilbert action. Under this generalized action, we attempt to explain 
i) the flat rotation curves observed in spiral galaxies, which is usually attributed 



Department of Physics, Tezukayama University 



to the existence of dark matter, and ii) the contradicting observations of uniform 
cosmic microwave background and non-uniform galaxy distributions against red- 
shift. For the former, we attain the flatness of velocities, although the magnitudes 
remain about half of the observations. For the latter, we obtain a solution with os- 
cillating Hubble parameter under uniform mass distributions. This solution leads 
to several peaks of galaxy number counts as a function of redshift with the first 
peak corresponding to the Great Wall. 
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1. Introduction 

One of the most straightforward generalizations of the Einstein-Hilbert action 
is the addition of a square terms of the scalar curvature and Ricci tensor. This 
action was introduced as counter terms to regulate ultraviolet divergences of the 
Einstein theory.' 11 It was also used to obtain a bounce universe, which avoids the 
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initial singularity of the big bang cosmology. 

In the present work we consider this action in order to overcome the difficulties 
encountered by the standard Einstein theory in explaining certain astrophysical 
and cosmological observations. The famous astrophysical observations that appear 
to be in conflict with our expectations from the Newton's theory of gravity is the 
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flat rotation curves of spiral galaxies. The rotation velocities are constant as a 
function of the distance from the center of galaxy, while a naive application of the 
1/r 2 force law implies a decline in the velocity function. This observation is usually 
accounted for by the incorporation of dark matter.' 41 In the present work we attempt 
to obtain flat rotation curves from the generalized Einstein action without relying 
upon dark matter. As for the cosmological observations, the recent technological 
developments in observations have yielded contradictory results : highly uniform 
cosmic microwave background W and non-uniform galaxy distributions in the scales 
of hundreds of Mpc.' 6 ' 7 ' In order to resolve this difficulty, Morikawa introduced a 
non-conformal scalar field in the action and obtained a theory with oscillating 
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Hubble constant. According to this theory, the mass distribution of the universe 
is uniform at any moment (hence, uniform cosmic microwave background), but the 
galaxiy distributions as a function of redshift become periodic (hence, non-uniform 
large scale structure of the universe) due to the oscillation in the expansion rate of 
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the universe. In the present work, we follow the same idea under the generalized 
Einstein action without introducing the scalar matter. 

The paper is organized as follows. In sect. 2 we discuss rotation curves of 
spiral galaxies under the generalized Einstein action. In sect. 3 we consider galaxy 
number counts under the generalized Einstein action. A concluding remark is given 
in sect. 4. 



In this section, we consider a modification of Newtonian theory of gravitation 
under the generalized Einstein action and apply this force law to the rotation curves 
of spiral galaxies. The generalized Einstein action is given by 



where G is the gravitational constant, R is the scalar curvature, R^ u is the Ricci 
tensor, A is the cosmological constants, and g is the negative determinant of the 
metric tensor. Throughout this paper we follow the conventions of Ref. 9 and, in 
particular, the speed of light is set equal to unity. The Einstein equation in this 
theory is then written as, 



2. Rotation Curves of Spiral Galaxies 



under the Generalized Einstein Action 





Riiv - T, R 9iiv - Ast^ + ci J M „ + c 2 K fil/ = 



SttGT, 




where and are defined by 





3 



= R,^ - HRfiu -^{UR + R a pR aP )g^ + 2R a < 3 R mv p , (2.4) 

with 

R,^ = V M V,i?, OR = 9 a(5 R-ap ■ (2.5) 

In order to obtain the modification to the Newton's 1/r 2 law under this theory, 
we consider the following static, weak field limit: 

£00 = -(1 + 20) , (2.6) 



g tJ =5 tj (l + 2ij) , (2.7) 

where and ip are functions of spatial coordinates only. Note that corresponds 
to the gravitational potential field. The Ricci tensor and scalar curvature are then 
given by 

R 00 = - A , (2.8) 
Rij^idf^ + df^ + A^ij , (2.9) 

i? = 2A0 + 4A^. (2.10) 

Substituting these into the time-time component of Eq. (2.2), we obtain the equa- 
tion for and ip: 

A[0 + 2ci A + (4ci + 2c 2 ) A -0] = 4vrGp , (2.11) 
where p is the mass density. The trace of Eq. (2.2), on the other hand , gives 
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another equation for and ip: 



[1 + (6ci + 2c 2 ) A] A (0 + 2^) = -4ttGp . 



(2.12) 



We have to solve (2.11) and (2.12) simultaneously. For the point source with 



p = M5 (3) (f) , 



(2-13) 



we obtain the following solution: 



[10,11] 



GM 



1 l e ~ mir 4 e~ m2r 
r 3 r 3 r 



(2.14) 



where 



V> = GM 



1 e" mir 
3 r 



2 e~ m2r 

3 r 



1 



6ci + 2c2 



C2 



(2.15) 



(2.16) 



Note that the first term in (2.14) is the usual Newton potential and that the second 
and third terms in (2.14) correspond to its corrections. The third term in (2.14) 
is, however, undesirable, since its existence does not yield the usual attractive 
Newton's force law in the limit r — > 0. We thus set 



c 2 = . 



(2.17) 



As for the constant c\, a positive c\ gives 



= -GM 



1 1/ cos ur sin ur 

- + o +2/ 

r 3 V r r 



)] 



(2.18) 
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where y is an arbitrary constant and 



= -S=» (2-19) 
V oci 



while a negative ci gives 



= _GM 



rl le 

- + 



— mir 



r 3 r 



(2.20) 



Both cases give the correct r — > limit of the observed gravitational constant Go 
with 

Go = ^ • (2.21) 

Eq. (2.20), however, implies more decline of force than the 1/r 2 law as a function of 
r, and this is opposite to what is expected from the rotation curves of galaxies. We 
thus adopt the solution of Eq. (2.18) and calculate rotational velocities of typical 
spiral galaxies. Taking into account the exponential mass distribution implied by 
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the surface photometry data, we approximate a spiral galaxy by a truncated disk 



with mass density 



m - \ : ^ * (Z) \ ^ ~ ^ (2-22) 
, (r > r g ) 



where we take 



r g = 4.3a [kpc] , a = 2, 3, 4, 5, 6 , (2.23) 

for the truncation radius.' 131 The value of co was also taken from Ref.ll. Namely, 
we impose the condition that the total mass is 2.2 x 1O 1O M for a = 2, where M 



is the solar mass. The free parameter fi and y in Eq. (2.18) were chosen so that 
the velocity curve becomes flat. Our choice is 
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[2.24) 



kpcJ 

V = 3 . (2.25) 

We thus obtain the square of the rotation velocity as a function of the distance r 
from the center of galaxy by 



2tt 



v 2 (r) 



—r f dr' f de'r'e- r '/ a -^-M\r-f"\) , (2.26) 
Air J J dr 



where 0o is the Green's function of Eq. (2.18) with M — 1, and we only consider 
r on the plane of the galaxy disk. 

In Fig. 1, we show the results of numerical integration of Eq. (2.26) for each 
value of a with r in the range of typical galaxy radius, i.e., 10 — 100 kpc. For a 
comparison, the result for a = 2 under usual Newtonian gravity is also shown by 
the broken curve in the figure. As is clear from the figure, the rotation velocities 
become flat for r ^ 50 kpc, while that of Newtonian gravity steadily decreases 
as r increases. However, the magnitudes of velocities are about half of what are 
observed for each value of a in Eq. (2.23) for r ^ 50 kpc. 



7 



3. Large Scale Structure of Universe 
under the Generalized Einstein Action 



In this section we study the cosmology under the generalized Einstein action. 
The action we consider is given by Eq. (2.1) with C2 = 0. Assuming the homoge- 
neous and isotropic universe, we employ the Robertson- Walker metric: 

ds 2 = -dt 2 + a 2 {t) ( - d ^ + r 2 dd 2 + r 2 sin 2 dd^j , (3.1) 

and the energy-momentum tensor T^ v of the form 

= pg^ + (p + p)XJ^\J v . (3.2) 

The time-time component and space-space component of the Einstein equation, 
Eq. (2.2), are then respectively written as 



H 2 + K — — — Qci(-H + K 2 + QHH 2 - 2KH 2 + 2HH) = — p , (3.3) 
3 3 



and 

-2H - 3H 2 - K + A + 6 Cl (2dfH + 12HH + 9H 2 + 18HH 2 - AKH 
-2KH 2 - K 2 ) = 8tiG P , 



(3.4) 



where 

H=±, (3.5) 
a 

and the dots stand for time derivatives. We now introduce a few dimensionless 



parameters as follows; 

r = H Q t , (3.6) 
^0 = - , (3.7) 

Pc 

k^-L, (3.8) 

A„ S ^, (3.9) 

where Hq is the Hubble parameter at the present time, and the critical density p c 
is defined by 

Pc - ^Hl . (3.10) 
The density parameter Qq have two components: 

n = n No + n Ro , (3.ii) 

where Qn q and f2# are respectively non-relativistic and relativistic components. 
We can then write 

3i?Q po a a 

where we have set the present value of the scale factor ao to unity: 

a = a(0) = l. (3.13) 

By (3.6) ~ (3.9), (3.12) and the energy-momentum conservation 

P = - 3 (l) (P + P)> ( 3 - 14 ) 
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we finally reduce (3.3) and (3.4) to 



1 fa 
VIcq [a* 



^Ro 
q ~r a 



1 f a 



2 V a 



a:- 
a' a" 



or \ or 



-) -^ + Ao 



a 



-) 2 -s(- 4 

a J \ a 



(3.15) 



and 



,1111 



12c 



respectively, where 



, +2— + - +^-3A 
or a V a / or 



a'a'" a' a" fa" \ (a 
4— =--12^- + 3 — +3 - 



0T 

a 2 I a 2 



a 



+ 4- 



(3.16) 



co 



cii/o 2 , 



(3.17) 



and the prime stands for derivative with respect to r. These are third-order 
and fourth-order differential equations for a(r). In particular, we consider solv- 
ing Eq. (3.16) numerically. For this we need initial condition of a,a', a" , and a'". 
We have chosen ao = 1 (Eq. (3.13)). It then follows that 



a = 1 , 



(3.18) 



a = -g 
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(3.19) 



where go is the deacceleration parameter 



go - -4? • (3.20) 

By choosing values of Qo,&o,Ao, and go which are observationally acceptable, we 
can determine a " for each value of cq from Eq. (3.15) evaluated at the present 
time: 

n - k + A + 6c [2a " - (a ') 2 + 2a ' — 3 + A;g — 2A; ] = 1 . (3.21) 



In the present work, we set 

k = , (3.22) 
A = , (3.23) 

Q Ro = 1(T 5 , (3.24) 

for simplicity and use only positive go. (Non-zero values of ko and Ao or negative 
go gave similar results unless their values are taken to be unnaturally large.) This 
leaves fijv , go(= ~ a o), an d co as free parameters. In Fig. 2 we show the results of 
numerical integration of Eq. (3.16) for several sets of values of these parameters. 
For co > 0, we have a de Sitter-like expansion (Fig. 2(a)), a bounce (Curve / in 
Fig. 2(6)) which is essentially the bounce solution of Ref. 2, an expansion followed 
by an eventual collapse (Fig. 2(c)), or an expansion, shrink, and eventual de Sitter- 
like expansion (Curve c in Fig. 2(d)). For co < 0, on the other hand, we generally 
have oscillation: oscillating expansion or oscillating bounce (Fig. 2(e) and (/)). 
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This last solution, Curve e in Fig. 2(e) in particular, is of interest, since it might 
explain the observed large scale structure of galaxy number counts as a function of 
red shift.' 6 ' 71 The number of galaxies dN which are located between the comoving 
distance r and r + dr at a fixed solid angle dfl with absolute luminosity between 
L and L + dL is given by 

dN = n(L, t)c?r 2 drdttdL , (3.25) 
where n(L, t) is the number density at time t with absolute luminosity L, satisfying 

n(L,t) = (^) 3 n(L,0) . (3.26) 
Since the redshift z and scale factor a are related by 

z = ^-l, (3.27) 



a 



we have 



°<*- MO + ltf. ,3.28) 



a J ao 

t t 



By substituting (3.26) and (3.28) into (3.25), we have 
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dN 

oc 



f t °(z(t') + l)dt'} 2 



(3.29) 



dzdQdL H(t) 
Note that Eq. (3.27) gives a differential equation for z(t) : 

^ = -(l + z)H(t) . (3.30) 

Hence, by solving Eqs.(3.16) and (3.30) simultaneously, we obtain the number 
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count in Eq. (3.29). The results for 



Q No = 0.1 , (3.31) 

a'o = -1.0 , (3.32) 

Co = -0.0002 , (3.33) 

are shown in Fig. 3 together with the scale factor a(r). Note that there exist 
four peaks in the galaxy count for z < 0.5 with the first one corresponding to 
the location of Great Wall! 61 Even though the peaks are not separated with equal 
interval, overall features are remarkably similar to observations. 



4. Concluding Remark 

In this paper we discussed a possible generalization of the Einstein theory by 
adding a square term, c\R 2 , of scalar curvature to the action. Choosing appropriate 
values of the coefficients, we found successful explanation of the flat rotation curves 
of spiral galaxies (c\Hq ps 5 x 10 -11 ) and the large scale structure of universe 
{c\Hq ps —2 x 10~ 4 ) at least at a qualitative level. Our theory, however, is 
not all faultless. It turned out that the coefficient 7 of the Robertson expansion is 
1/2, and this disagrees strongly with, for instance, the radar echo experiments from 
Mercury, which implies that 7 is very close to 1 in accord with the Einstein theory! 141 
Originally R 2 term was introduced as a result of quantum corrections, 111 so that the 
coefficient c\ will have the logarithmic r dependence as c\ = (const.) x log(r/ro). 
This interpretation may lead to a consistent explanation of both the rotation curves 
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of spiral galaxies and the large scale structure of the universe from the single 
effective Lagrangian with the running coefficient c\ that depends on the distance 
scale of interest. The difficulty of 7 = 1/2 may also be avoided if tq takes a value 
of the order of solar radius. 

Acknowledgements: We are grateful to M. Fukugita, H. Kawai, K. Kume, and K. 
Uehara for useful discussions. 
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FIGURE CAPTIONS 

1) Rotation curves of spiral galaxies under the generalized Einstein action. The 
numbers on the curves indicate the values of a in Eq. (2.23). The prediction 
of Newtonian theory for a = 2 is also given by a broken curve for comparison. 

2) Scale factor of universe as a function of r for various values of parameters 
co, f2 No j and a ' . The parameter cq is positive for (a) ~ (d) and negative for 
(e) and (f). The values of |co| are indicated in the figure by the letters a, b, 
c, d, e, / for 10°, 10 _1 , KT 2 , 1CT 3 , 1CT 4 , 10 -5 , respectively. 

3) (a) Scale factor of universe as a function of r for the set of parameters in 
(3.31) ~ (3.33) and (b) number count of galaxies as a function of z for the 
same set of parameters. The scale of the ordinate for (b) is arbitrary. 
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